Non-adiabatic phonons within the doped graphene layers of XC6 compounds 
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We report Raman scattering measurements of BaC6, SrC6, YbC6, and CaC6, which permit a 
systematic study of the phonons and the electron-phonon interaction within the doped graphene 
layers of these compounds. The out-of-plane carbon phonon softens as the spacing of the graphene 
layers is reduced in the series BaC6, SrC6, YbC6, and CaC6- This is due to increasing charge in 
the 7r* electronic band. Electron-phonon interaction effects between the in-plane carbon modes 
at ~ 1500 cm -1 and the 7r* electrons cause a strong non-adiabatic renormalization. As charge is 
transferred into the 7r* band these non-adiabatic effects are found to increase concurrent with a 
reduction in the phonon lifetime. 

PACS numbers: 71.20.Tx,63.20.kd,74.25.Kc,78.30.-j 



I. INTRODUCTION 

There is mounting evidence that the nature of electron- 
phonon interaction in graphitic systems is not fully 
understood (e.g. in Graphene^ - — or carbon nanotubes 
(CNTs) 4 ). In graphite electron-phonon scattering is 
the main contribution to the phonon lifetime and there- 
fore the Raman linewidth, which is typically 11.5 cm -1 . 5 
When atoms are inserted into graphite, electrons often 
act to dope the graphene sheets and the Raman linewidth 
increases dramatically to ~ 100 cm -1 . The question 
of the electron-phonon interaction is also interesting in 
light of the discovery of superconductivity at 11.5 K in 
CaC6i&£ Measurements of a finite Ca isotope effect^ and 
of an s-wave symmetry superconducting gap^^ suggests 
that the superconductivity is electron-phonon mediated. 
However, several experimental and theoretical techniques 
implicate different groups of phonons and electrons are 
responsible for the superconductivity^! 1 - 1 - - — 

The adiabatic Born-Oppenheimer approximation, in 
which electrons are assumed to instantaneously respond 
to the motion of the ions, has been widely applied within 
the framework of density functional theory (DFT) to 
predict the phonon frequencies in metals^ Despite this 
approximation not being theoretically justified, the dis- 
crepancies introduced are usually small < 1 8 How- 
ever, a much larger effect can be observed provided 
\q-v F \ <C o V 8 ' 19 where q is the phonon wavevector, 
vp is the Fermi velocity, and cj is the phonon fre- 
quency. Graphene and materials composed of graphene 
layers typically have small Vp perpendicular to the 
graphene sheets. This makes these systems ideal to ob- 
serve non-adiabatic effects. Weak non-adiabatic effects 
have already been observed in metals such as osmium, 20 
graphene,— graphite,— and CNTsj^ Graphite interca- 
lation compounds (GICs) provide an ideal opportunity 
to examine phonons within the heavily electron-doped 
graphene sheets in these compounds, where the increased 
electronic density of states at the Fermi surface should 



TABLE I: The values of d and the T c for the GICs under 
consideration.—^ ' 14 ' 23 





d(A) 


T, (K) 


BaC 6 


5.25 


<0.08 


SrC 6 


4.95 


1.65 


YbC 6 


4.57 


6.5 


CaC 6 


4.52 


11.5 



make these effects more significant. 

Recent theoretical work has shown that for most GICs 
the completely non-adiabatic zone-center phonon fre- 
quencies oj na are closer to the experimental values than 
the adiabatic frequencies ui A M- Using u A and u> NA along 
with the experimental phonon frequency ujbwf, Saitta 
et al. calculate the contribution of the electron-phonon 
scattering to the phonon full width at half maximum 
7^ p , which can be compared to the experimental width 
r. This formalism works quite well for KCs and CaCg, 
but the width is overestimated in RbCg. Furthermore 
LiCg and KC24 fall outside the boundaries of uj a and 
lu na . To the best of our knowledge, we present the first 
systematic experimental study of these effects within a 
particular family of GICs and the first published phonon 
spectra of BaC6, SrC6 and YbC6. The solid experimen- 
tal trends presented here, and the interpretation thereof, 
can be applied widely within graphitic materials as the 
phonons we are measuring are the direct relative of the G 
mode in graphite, graphene, and CNTs. We choose sam- 
ples with the same in-plane superlattice: BaCg, SrCg, 
YbCg, and CaC6- Changing the intercalant within these 
GICs controls the separation of the graphene sheets d, 
which provides a clear physical parameter linked to the 
changes in T c in these compounds. Table fl] shows that as 
d decreases the superconducting transition temperature 
T c increases dramatically. We show that reducing d, by 
changing the intercalant, leads to a softening of the C 2 
phonons modes, which is attributed to charge transfer 
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FIG. 1: (Color online) The Raman spectra from different po- 
sitions numbered on the GIC surface shown in the increased- 
contrast inset. Spot 1 corresponds to pure, ordered GIC with 

1500 cm" 1 . Spots 2-4 also show 
1600 cm" 1 phonons which arise 
from sample disorder. All spectra are normalized to represent 
the same counting time. 



one C xy mode visible at : 
the ~ 1350 cm" 1 and the 



FIG. 2: (Color online) The Raman spectra of the C z modes for 
BaC6, SrC6, YbC6, and CaC6- Black dots represent the data 
points and the solid lines are the Lorentzian fits. A smooth 
background has been removed from the data. The inset shows 
the variation in peak position for the experimental frequency 
u>o (dots) and the DFT value lj a (crosses)^. 



into the ir* states. Using this, we go on to demonstrate 
that the C xy phonon modes are strongly non-adiabatic, 
and that the size of these effects increases with increas- 
ing charge transfer to the n* states. Furthermore, our 
work provides new evidence that the size of the electron- 
phonon interaction is stronger than the current theory 
predicts^ 



II. EXPERIMENTAL METHODS 

Samples were synthesized using the vapor transport 
method&i^ from natural Madagascan graphite flake. X- 
ray diffraction confirmed the phase purity of the samples, 
which were typically ~ 200 fim in size. For compari- 
son CaC6 and BaCg were also synthesized from highly 
oriented pyrolytic graphite (HOPG) using the Li alloy 
method.— The spectra were measured at room tempera- 
ture using a Renishaw In Via confocal micro-Raman sys- 
tem, equipped with a 514.5 nm Argon ion laser, which 
was focused to ~ 3 [im onto the ab-plane of the sample. 
The laser power at the sample was kept below 4mW to 
avoid laser heating effects. 

When dopant atoms are introduced into graphite, the 
in-plane C xy mode of E2 S symmetry (also called the G 
mode), originally at ~ 1580cm" 1 in graphite and CNT, 
softens and increases in linewidth. The new superlattice 
also folds a mode present at the K-point in the graphite 
Brillouin zone to the T-point, making it Raman active. 
This C z mode, around 500 cm" 1 , involves out-of-plane C 



motion. 

Alkali and alkali-earth GICs react readily with oxygen 
and wateri 25 i 26 The samples were therefore cleaved in an 
argon glove-box (H2O and O2 levels <0.1ppm) before 
being immediately transferred to pre-baked quartz tubes 
and sealed under a vacuum of < 10 -6 mbar. 

The inset of Fig. [T] shows a typical increased contrast 
image of the surface of the GICs considered here. While 
to the naked eye all samples look extremely shiny, under 
high magnification one can see and avoid step-edges or 
small areas of surface degradation, which cause the sur- 
face to appear slightly blackened. This inhomogeneity 
explains why Refs. 1271)281 observe contaminant peaks in 
most of their CaCe spectra. Only by choosing the most 
reflecting and therefore the most pure part of the sample 
(position 1 in Fig. [TJ is a clean spectrum obtained con- 
taining a C xy mode at w 1500 cm" 1 and a C z mode at 
« 450 cm -1 (too weak to see in this plot). We assign, 
in agreement with Ref. 127J the peaks at ps 1350 and ~ 
1600 cm" 1 to disordered graphite to which small areas of 
the surface degrade. These contaminant peaks are inher- 
ently much more intense than the true GIC Raman peak, 
because Raman scattering from disordered graphite is al- 
ways a resonant process i 29 ' 30 Ref [H attributes the peak 
at ?» 1350 cm" 1 to disordered CaC6- We note that this 
mode is not associated with the pristine surface. 
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FIG. 3: (Color online) The Raman spectra of the C xy modes 
for CaC6, YbC6, SrCe, and BaC6- Black dots are the mea- 
sured spectrum, the thin line is the BWF function, the dotted 
lines are the surface contaminant and the fit is shown as a solid 
line. A linear background has been removed from the data 
and the spectra have been offset from one another for clarity. 



III. RESULTS 

The C z modes from BaC6, SrC6, YbC6, and CaCg are 
shown in Fig. [5J A clear trend in the peak positions 
is observed, with the peak softening from 480(3) cm -1 
in BaC 6 , 470(3) cm- 1 in SrC 6 , 443(3) cm" 1 in YbC 6 to 
438(3) cm -1 in CaC6. To date, experimental phonon 
measurements of this family of GICs has been limited 
to CaCg. Out-of-plane phonon dispersions below 50meV 
have been measured by inelastic x-ray scattering^!, and 
found to be close to the DFT predictions. Higher en- 
ergy modes above 50meV have only been measured 
via Raman scattering, for CaCg Hlinka et al. report « 
450 cm - X 2L and Mialitsin et al. obtain 440 cm - 1 

The C xy modes in GICs often possess the asymmetric 
Breit-Wigner-Fano (BWF) line shape^i^i 3 ^ 3 - which is 
modeled as a Raman signal intensity of 



I 1 T qT/2 > 

1 + ( u "r/T F ) 2 



(1) 



where 1/q quantifies the asymmetry of the shape and 
wbwf and T are fitting parameters to the central fre- 
quency and full width at half maximum respectively. In 
the limit 1/q — > a Lorcntzian line shape is recovered 
with a full width T and frequency luq. Typical Raman 
spectra of the C xy modes for the samples considered here 
are shown in Fig. [3] The BWF line shape (Eq. [T]) was 
used to fit the peak and Lorentzian functions were also 



TABLE II: The fit parameters for the C xy mode of the GICs 
considered here. Values in columns ui A and w NA correspond 
to the DFT calculations of Ref. 18i All values except 1/q are 
in cm" . 





OJBWF 


r 


1/q iv A uj NA a 




BaC 6 


1510(5) 


70(10) 


-0.25(5) 1462 1521 700(200) 


46(8) 


SrC 6 


1510(5) 


100(10) 


-0.30(5) 1459 1530 910(160) 


64(5) 


YbC 6 


1515(5) 


100(10) 


-0.30(5) 




CaC 6 


1510(3) 


120(10) 


-0.35(5) 1446 1529 790(80) 


70(4) 



included in some fits to account for any small contamina- 
tion from the D and G modes. The errors were estimated 
by fitting a number of different spectra. Table |TT] pro- 
vides a list of the fitting parameters. The data for CaC6 
are in excellent agreement with Hlinka et al. who report 
0^^=1511 cm" 1 , r^lllcm- 1 and 1/q = 0.28,22 Mi- 
alitsin et al. measure a spectrum at 20 K showing large 
D and G peaks^ They describe a Lorentzian line shape 
with wo=1508cm _1 and r=76cm~ 1 . 



IV. ANALYSIS AND DISCUSSION OF RESULTS 

The trend in the C z modes to soften from BaCe, SrCg, 
YbCg, to CaCg clearly depends on d rather than the in- 
tercalant mass. We propose that by increasing d charge 
is transferred out of the 7r* band depopulating the Dirac 
cones. As the ir* band is anti-bonding filling it with elec- 
trons destabilizes the bonds causing the softening of the 
phonons, thus these modes are sensitive to charge trans- 
fer. Boeri et al. predict C z modes soften considerably 
when charge is transferred into the it* states in jellium 
intercalated graphite^ The inset of Fig. [5] compares the 
peak values to adiabatic DFT calculations. 13 Theory pre- 
dicts that the difference between u> A and u NA is due 
to intra-band electron-phonon scattering, as inter-band 



process contribute equally to u> A and iu NA & Given that 
the experimental values are close to uj a , and assuming 
that the Ig.f^l <C w condition is fulfilled, this agree- 
ment shows that there is very little renormalization of 
these modes due to intra-band transitions. In addition 
the line-widths of these modes are relatively small and, 
within the experimental errors, show no trends with d. 
The strong doping sensitivity of the C 2 modes means that 
the experimental frequency can be used, in conjunction 
with DFT, to infer a bulk value of the charge transfer. 
The experimental value for CaCg of 438(3)cm _1 is close 
to the DFT prediction of 433 cm -1 , 13 In this model 0.20 
electron per C atom2£ are found in the ir* states, which 
is close to the ARPES derived value of 0.22J£ This sup- 
ports the view that the charge transfer in the surface lay- 
ers probed in angle-resolved photoemission spectroscopy 
(ARPES) is representative of the bulk, and shows that 
these techniques are consistent in this regard. The charge 
transfer dependence proposed here also explains the se- 
ries of GICs CsCs, RbCs to KCsr 5 which show a similar 



4 



trend. The gradient of these effects is similar, but the 
overall values are different, since in XCg GICs (where X 
= intercalant) the charge transfer is less. We also note 
that in XCs GICs the C z mode is folded from M in the 
graphite Brillouin zone, whereas for IC6 GICs the mode 
is folded from K, which is at slightly lower energy. 

The uibwf values shown in table IPTl show no depen- 
dence on intercalant mass: these phonons only sample 
the environment of the (doped) graphene planes. These 
values are all significantly different to the values pre- 
dicted using DFT— with the adiabatic approximation 
uj a (see Table HI]) . In fact, the values are closer to the 
completely non-adiabatic limit oj na , which demonstrates 
that non-adiabatic effects are essential to understand 
these phonon energies. The degree of non-adiabaticity, 
may be estimated by u>bwf — uj a which increases as 
charge is transferred to the ir* states in BaCg, SrCg to 
CaCg. Despite the similarity in ujbwf, T varies signif- 
icantly increasing in the order BaCg, SrCg f» YbCg to 
CaCg. This is indicative of a decrease in phonon lifetime 
as charge is transferred into the it* states. The trends 
observed when doping the it* states by changing the in- 
tercalant atom in GICs are different to those measured 
in graphene and metallic CNTs at the much lower dop- 
ing levels accessed by electrostatic gating. 1 These sys- 
tems show an increase in phonon energy and a decrease 
in r with electron-doping. However, different effects arc 
present in the high doping case of the GICs because the 
Kohn anomaly which exists near q = in graphene 3 - is 
displaced to finite q and the tt* band is altered and is ac- 
companied by other intercalant-related bandsi 11 ' 13 ' 25 ' 34 

Within Saitta et al.'s theory ujbwf should fall in be- 
tween the completely adiabatic limit uj a and the com- 
pletely non adiabatic limit uj na . The electron scatter- 
ing rate a = H/t, where r is the electron momentum- 
relaxation time controls were ujbwf falls in between 
these limits. For a 3> u> A we find ojbwf ~> u> A and 
for a <C oj a the system becomes fully non-adiabatic and 
ujbwf u> NA . We observe no systematic trend in a, 
and we find the values are highly sensitive to any small 
errors in the experimental and theoretical frequencies. 
It is therefore difficult to use this method to accurately 
extract the electron scattering rate. It is however, pos- 
sible to rearrange Saitta et al.'s equations for derive the 
electron-phonon scattering induced linewidth 



— ~ y (ujbwf - uj a )(uj na - ujbwf)- (2) 

The derived values of ^ PC account for the trend in T 
to increase going from BaCg, SrCg to CaCg, but sig- 
nificantly underestimate its absolute value. The F val- 
ues reported here should be the intrinsic linewidth in 
these GICs, with a negligible disorder contribution, as 
these results are reproducible in samples produced from 
Madagascan flake graphite and HOPG. This is despite 
HOPG having a crystallite size of ~ 1 /im, far smaller 
than the ~ 100 [im in Madagascan flake graphite. An 
order of magnitude calculation can be performed to esti- 



mate the length scales which correspond to the observed 
r. The C xy modes in these GICs have a typical lifetime 
of ujbwf /r w 15 periods and stretch bonds over ~ 5 A, 
which implies the atomic motion is correlated over only 
- 80 A. Mada gascan flake based CaCg and YbCg sam- 
ples made in the same way as the samples used in this 
study have ordered regions of ~ 50 Furthermore, 
the mean distance between scattering centers, as esti- 
mated from the low temperature resistivity is ~ 1000 A.— 
Disorder broadening should therefore have an insignifi- 
cant role here. Anharmonic effects are also unlikely to 
cause such a large linewidthJ^ We therefore suggest that 
theory underestimates the degree of electron-phonon in- 
teraction in these compounds^ 

These results also have implications for the debate 
outlined in the introduction (Sec. [J) regarding which 
phonons and electrons are most relevant for the super- 
conductivity in SrCg, YbCg and CaCg. The trend un- 
covered that T c increases when charge is transferred into 
the graphene layers provides strong evidence that the C 
phonons and the tt* electrons are involved in the super- 
conductivity. This is very difficult to reconcile with the 
isotope effect measurements^ and early theory^ which 
suggest that superconductivity is predominately due to 
the intercalant atoms. Our work supports the assertion 
that the C atoms are relevant, either the C z mode o 11 ' 13 ^ 4 
or the C xv modest 



V. CONCLUSIONS 

In conclusion, we present phonon measurements of 
BaC 6 , SrC 6 , YbC 6 , and CaC 6 . Raman spectroscopy 
is employed to monitor the systematic changes in both 
the phonons and the electron-phonon interaction in the 
graphene sheets as d is varied by changing the inter- 
calant atom. The C z phonons are found to soften in the 
order BaC6, SrCg, YbCg to CaC6 indicative of change 
transfer into the tt* states as d is reduced. Concurrent 
with this charge transfer is a large increase in T c . While 
the out-of-plane carbon mode is well described within 
the Born-Oppenheimer approximation, the in-plane car- 
bon mode is not. We demonstrate that the in-plane 
carbon modes show strong non-adiabatic effects, which 
lead to phonon frequencies that are considerably higher 
than the adiabatic prediction. As the charge transfer to 
the tt* states increases, these effects get stronger caus- 
ing an increasing discrepancy between the adiabatic pre- 
diction and experimental frequency and an increase in 
the phonon linewidth. Theoretical work (Ref. [l8l ) an- 
ticipates the relevance of the non adiabatic effects in 
these GICs. However, our work provides new evidence 
that the degree of electron-phonon interaction for this 
mode is stronger than the current theory predicts. Non- 
adiabatic effects can be used to extract information on 
the electrons, phonons and the electron-phonon interac- 
tion within the doped graphene layers, which constitute 
these compounds. Through a systematic study of XC6 



5 



GICs, we identify trends which should be widely applica- 
ble to characterize a range of graphitic compounds using 
Raman spectroscopy. 
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